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Abstract 

The spin network quantum simulator relies on the su(2) representation ring (or its 
g-deformed counterpart at q — root of unity) and its basic features naturally in- 
clude (multipartite) entanglement and braiding. In particular, g-deformed spin net- 
work automata are able to perform efficiently approximate calculations of topolog- 
ical invarians of knots and 3-manifolds. The same algebraic background is shared 
by 2D lattice models supporting topological phases of matter that have recently 
gained much interest in condensed matter physics. These developments are moti- 
vated by the possibility to store quantum information fault-tolerantly in a physical 
system supporting fractional statistics since a part of the associated Hilbert space 
is insensitive to local perturbations. Most of currently addressed approaches are 
framed within a 'double' quantum Chern-Simons field theory, whose quantum 
amplitudes represent evolution histories of local lattice degrees of freedom. 

We propose here a novel combinatorial approach based on 'state sum' mod- 
els of the Turaev-Viro type associated with S(7(2) ? -colored triangulations of the 
ambient 3-manifolds. We argue that boundary 2D lattice models (as well as ob- 
servables in the form of colored graphs satisfying braiding relations) could be con- 
sistently addressed. This is supported by the proof that the Hamiltonian of the 
Levin-Wen condensed string net model in a surface E coincides with the corre- 
sponding Turaev-Viro amplitude on £ x [0,1] presented in the last section. 
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1 Braiding and entanglement in spin networks 



By spin networks we mean here 'computational' graphs the nodes and edges of which 
are labeled by dimensions of SU (2) irreducible representations (irreps) and by SU (2) 
recoupling transformations, respectively. For this reason spin networks can be thought 
of as a generalized quantum computational framework for dealing with unitary trans- 
formations among (entangled) many-angular momenta states. Actually, the compu- 
tational space of the spin network simulator^ encodes in its very definition the rep- 
resentation ring of the Lie group SU (2) given by finite-dimensional Hilbert spaces 
supporting irreps of SU(2) endowed with two operations, tensor product <g> and direct 
sum © (providing a ring structure over the field C) as well as unitary operators relating 
(multiple tensor products of) such spaces. 

Unlike the usual standard circuit model[2|, here it is possible to handle eigen- 
states of N (pairwise coupled) angular momenta labeled by integers and half-integers 
ji , j2 , • ■ ■ ! Jn (in h units) and not simply 7V-qubit states for which j\ = j2 = ■ ■ ■ = 
jjf = 1/2. The (re)coupling theory of N SU(2) angular momenta provides the whole 
class of unitary transformations that can be performed on many-body quantum sys- 
tems described by this kind of states[3|. Any such general transformation -expressed 
in terms of a 3nj-coefficient (iV = n+ 1) when a choice of basis sets is made explicit- 
can be split into a finite sequence of two basic unitaries, called twist (or 'trivial braid- 
ing') and associator in the language of tensor categories^]. The former acts on the 
(ordered) tensor product of two Hilbert spaces V, W supporting irreps of SU (2) by 
swapping them, namely 

Tvw ■ V <g> W -> W <E> V with T wy oT V:W = ld v ® w (1) 

and its action on a binary coupled state amounts to a trivial phase transform. 
The associator F relates different binary bracketing structures in the triple tensor prod- 
uct of irreps V, U, W (intertwiner) 

F : (V ® U) <g> W -> V ® (U <g> W) (2) 

and is implemented on any binary coupled state as a transformation involving one 
Racah-Wigner 6j symbol. Notice that both ([TJ and are isomorphisms but the 
associator reflects a physically measurable modification of the way in which inter- 
twiner spaces are coupled. More complicated multiple tensor products can be related 
by various combinations of braidings and associators, so that the basic isomorphisms 
must satisfy compatibility conditions, a so-called pentagon identity and two hexagon 
identies[4| (in SU(2) angular momentum theory they correspond to the Biedenharn- 
Elliott identity and the Racah identity, respectively[3|). 

The most effective way of dealing with non-trivial braiding operators -to be used in 
connection with the study of both braid group representations and braid statistics (typ- 
ically occuring in anyonic systems)- is to move to the representation ring D\ (SU{2) q ) 
(modular tensor category) of the q-deformed Hopf algebra of the Lie group SU{2), 
SU{2) q (q a root of unity). Then the operation T in ([TJ turns out to be substituted by a 
non-trivial braiding 

Rv.w ■ y ® W -> W <g> V with R wy o R VM ^ 1& V ® W . (3) 

The q- deformed spin network model of computation! 5 1 is modelled on the ^-representation 
ring 

(X(SU(2) q );R;F), (4) 
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where F denotes from now on the ^-counterpart of the associator (f2]). Once suitable 
basis sets are chosen, the isomorphisms R and F can be made explicit. In particular, 
F contains the g-deformed counterpart of the 6j-symbol and, regarding it as a uni- 
tary matrix, it is also referred to as 'duality' (or 'fusion') matrix in conformal field 
theories! 6j . 

The framework outlined above was exploited in a series of papersE 13, where 
families of g-deformed spin network automata were implemented for processing effi- 
ciently classes of computationally-hard problems in geometric topology -in particular, 
approximate calculations of topological invarians of links[8 | (collections of knots) and 
3-manifold[9|. A prominent role was played there by 'universal' unitary braiding oper- 
ators associated with representations of the braid group in D\ (SU (2) g ). Suitable traces 
of matrices of these representations provide polynomial invariants of 5t/(2) g -colored 
links, while weighted sums of the latter give topological invariants of 3-manifolds 
presented as complements of knots in the 3-sphere. These invariants are in turn rec- 
ognized as partition functions and vacuum expectation values of physical observables 
in 3-dimensional Chern-Simons-Witten (CSW) Topological Quantum Field Theory 
(TQFT)[10|. The CSW environment actually provides not only a physical interpreta- 
tion of such quantities but it is universal in the sense that includes, besides the quantum 
group interpretation quoted above, monodromy representations of the braid group aris- 
ing in a variety of (boundary) conformal field theories [6| (where point-like excitations 
confined in 2-dimensional regions evolve along braided worldlines). 

It is worth mentioning that the g-spin network approach represents a naturally 'dis- 
cretized' version of the topological framework for quantum computation proposed a 
few years ago and further improved recently! 1 1 1- in the light of some basic questions 
raised by that paper, we are going to outline a novel approach to the whole matter 
of topogical phases (section 2) providing in particular a transparent combinatorial de- 
scription of condensed strings nets (section 3). 

2 <S'C/(2) g -colored triangulations: 3D invariant parti- 
tion functions and induced 2D lattice models 

The issue of connections between (topological) gauge theories in 3 spacetime dimen- 
sions and 2D (integrable) lattice models! 12 | nas been intensively investigated over the 
years in a variety of different contexts. The renewed interest driven by the search for 
a fault-tolerant quantum computer based on manipulations of Non-Abelian quantum 
Hall states | TT][T3) represents a major challenge on both experimental and theoretical 
grounds. Most of the currently addressed approaches are framed within a CSW envi- 
ronment, but an ab initio discretized framework could reveal much more effective, as 
briefly outlined in the following. 

The Turaev-Viro (TV) 'state sum' model lfT4l provides a well-defined {i.e. finite) 
topological invariant for any closed 3-manifolds M 3 . Given a triangulation T 3 of A/ 3 , 
namely a dissection into tetraheda, an 5[/(2) 9 -coloring is assigned to it according to a 
set of admissible initial data (not made explicit here). A state functional for T 3 is built, 
where the Ao vertices, Ai edges and A3 tetrahedra are suitably weighted. Finally, a 
summation over all admissible colored triangulations of M 3 is performed, the resulting 
functional being invariant under a set of combinatorial moves ensuring that its value 



3 



depends only on the topological type of the manifold M 3 . The TV state sum reads 

N ± N 3 ( . . . . 

{j} A=l B=l ^ 4 J5 J6) B 

where the summation is over all colourings {j} labeling irreps of SU (2) q 
(q = exp{27ri/r}, with {j = 0, 1/2, 1 . . . , r-l»; w A = (-l) 2 ^ [2 ]A + l] q where [], 
denote the quantum dimension of the irrep; w = 2r j (q — q^ 1 ) 2 and { :::}b represents 
the q — 6j symbol whose entries are associated with the six edges of tetrahedron B. 

The invariant © equals the square modulus of the Reshetikhin-Turaev invariant|9|, 
which in turn represents the Chern-Simons partition function Z cs f° r an oriented 3- 
manifold M 3 , namely Z T v [M 3 ; q] < — ► | Z C s [M 3 ;k] | 2 , where k — 2(r — 1) is 
the level of Chern-Simons functional. This feature is crucial in view of applications to 
topological phases, where a 'doubling' of the basic CS setting is needed lfTTI . 

The TV model can be further extended to deal with both 3-manifolds with 2D 
boundary component(s) and 'observables' related to embedded links and (ribbon) graphs 
The extension to a pair (M 3 , E) with E = dM 3 can be carried out in two different 
ways, depending on whether the boundary surface E inherits a fixed triangulation lfT5l 
or a 'fluctuating' one[ 16 |. Observables in the form of colored graphs, satisfying braid- 
ing relations, can be consistently introduced lfl5l [TTIl in any oriented triangulated com- 
pact manifold (A/ 3 , E). 

On the basis of such results we are currently addressing a combinatorial reformula- 
tion of theoretical foundations of topological phases lfl8l . A first result is discussed in 
the last section. 

3 Topological spin liquids: a combinatorial description 

We will now briefly introduce the condensed string nets of Levin and Wen lT9ll in (2+1) 
dimension and show that the projector to the ground state of the boundary model can 
be obtained from the corresponding Turaev-Viro partition function Ztv(E x [0, 1], X) 
for the cylinder M 3 — E x [0, 1] (E is oriented and X denotes a fixed identical trian- 
gulation on both E x {0} and E x {1}). 

Consider the string net model defined on the honeycomb lattice lfl3ll with micro- 
scopic degrees of freedom associated to the edges. The degrees of freedom are called 
string types and are endowed with N + 1 labels and orientations. The change of ori- 
entation of an edge is equivalent to the change of the label i to its dual i* and the 
label is self-dual. The Hamiltonian is a sum of potential and kinetic energy. When the 
relative coefficient of the kinetic piece is large, the ground state is a condensate of a 
dense string net. The universal long distance behaviour of these nets is characterized 
by the so-called fixed-point wave functions on the space of all configurations. The 
fixed-point wave functions $ are in one-to-one correspondence with modular tensor 
categories ll20ll and obey the following rules (written by using the original notation |fl9l . 
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although slightly different from the conventions in (|5]l): 



<T> | j =* -Ljj- j (6) 

$ ( O') =rfi$ C ) (7) 



<!>( A I =,),,<!> I T C> T | ,8) 



1 i j ; _u y 

* ( - e (KG!) 

where the F matrix and the quantum dimensions di satisfy consistency conditions 
given by 



F ijk _ J!*_x. 



u 3 

pijm plkm" pjim pimj v mVn 

kin r jin Ikn* k'nl 



VjVl 
N 

Ejpmlq pjip pjs'n _ pjip priq* 
kp*n mns* r lkr* q* kr* mis* 

n=0 

with Vi = \fd~i (= w,; in TV notation) and Sijk = 0, 1. The latter is the branching rule: 
the triple {ijk} is admissible at a vertex when S^k = 1. These rules coincide with 
those of Turaev and Viro lfT4l for the evaluation of closed manifold invariants in three 
dimensions and the (first two lines of the) consistency conditions encode the standard 
properties of the q — 6j symbol suitable normalized as 

I (ID 

while the last relation in (TTOt is the Biedenharn-Elliott (pentagon) identiy (cfr. section 
1). 

The Hamiltonian for the string-net model on the honeycomb lattice reads 

i P 

where / runs over vertices and p over plaquettes. Qj — S^k with i,j, k being the 
colours of the edges incident to the vertex /. The 'magnetic constraints' are given by 
the sum X^fLo a sBp with to be described below. We consider the case when the 
coefficients are 

(13) 



This normalization prescription corresponds to the existence of smooth continuum limit 
of the lattice model in one hand and gives the precise weight of the Turaev- Viro eval- 
uation on the other, as proved below. 
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Figure 1 : a) shows a part of the honeycomb lattice and its dual graph. In b) the 'plane' 
of the triangle abc is the 'plane' of the honeycomb lattice and there is another one 
with triangle a" b" c" . The matrix element of the projector is interpreted as a three 
dimensional extrapolation between the bra and the ket string net state given by a trian- 
gulation. The prism in the figure correspond to a part of the projector P where there is 
a tetrahedron for every associated F symbol. 



The magnetic constraints 



There is a nice intuitive way of determining the matrix elements of the magnetic con- 
straints described in appendix C lfl9ll . In short, for each hexagon of the lattice, one 
draws a smaller concentric hexagon with labels s at every edge in hexagon p, then (i) 
connects the hexagons p and s with new edges and labels them with (this does not 
change and (ii) use the rules (O repeatedly until one reaches the single hexagon 
again. The resulting coefficients define the matrix elements of the constraint given 
esplicitly by 





(14) 



Epbg'h pch'i pdi* j pej* k pfk'l j?al* g 
s*h' g'* s*V h" r s* j'i" s* k' j" s*l'k'* s*g'l'* 

g' h' i' j' k' I' 



The operators B p with the coefficients given by dT3t are projectors (as the Q/'s) and 
all these constraint operators commute with each other. 

The image of the projector P = Q q B p is the ground state of the model. 
We are going to show explicitly that 



P = Ztv(Z x [0,1], X) 



(15) 



where E is the (oriented) surface on which the honeycomb lattice is defined and X is 
the dual triangulation of S x {0} ~ E x {1} ~ E. To build the dual graph one places 
a vertex in the center of each face and connects the new vertices in neighbouring faces. 
In particular, the dual of the trivalent honeycomb lattice is a triangular graph as shown 
in Fig.[TJi). 

Let us start constructing J\ B p . Each operator B* contains six F symbols corre- 
sponding to the six vertices of a hexagon. We collect the three F symbols from three 
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different operators corresponding to a given vertex. Their order does not matter 
since the B operators commute. Choosing the order B^B^B^ (cfr. Fig.QJ) we get 
the contribution Ff c A ,F? a J. Ff^',- a '. Note also that we have omitted denoting the 

S\0 c S2C a 53a o o 

orientation for the sake of simplicity, but it is not difficult to incorporate it consistently. 
Comparing this expression with Fig. QJ)), we find that it corresponds to the triangula- 
tion of the prism with boundary triangles abc and a" b" c" . In this identification we have 
associated a q — 6j or F symbol to every tetrahedron, which is the prescription of the 
Turaev-Viro state sum 

The second step consists in checking whether the gluing along the triangulated 
boundary quadrilaterals can be made consistent with the algebraic pattern. Using the 
freedom of multiplying the F symbols at a given vertex in arbitrary order (due to 
the commuting B p operators that they are part of), the consistency can be achieved. 
It is clear now that we have a triangulation given by the dual graph with edge la- 
bels a" b" c" d" . . . and an identical one with labels abed . . . and these bound a 3- 
dimensional triangulation given by triangular prisms suitably glued together. 

We now check whether the weights associated to edges and vertices match that 
of the TV prescription. Taking into account the normalisation (fTTT i we find that there 
is a factor of di associated to each internal edge denoted by a'b'c'd 1 .... Because of 
the coefficient ( fT3l l it is also true for the internal edges 'perpendicular' to E labelled 
by Si in the figures (this weight comes from the numerator of a s , the denominator 
will be associated to vertices bounding these edges). To the boundary edges with la- 
bels a" b" c" d" . . . there is also a d a " , (4» etc. associated, whereas to those with label 
abed . . . there is no non-trivial weight. This is consistent with the TV prescriptiorQ. 
The last type of simplexes we have and not discussed yet are the boundary vertices. 
We associate the square root of the denominator Q]\ d 2 )~ x / 2 = D^ 1 (= w _1 in TV 
notation) of the coefficient a s to the two boundary vertices connected by the edge la- 
beled by s. This way we covered both boundaries and found agreement with the TV 
prescription. 

The last step is the sum over internal colouring, which has to be performed to 
get the partition function. This corresponds to the summation over the labels Sj and 
a'b'c'd' . . . , which come from ^ a s B p and the matrix multiplication in J\ B p . Fi- 
nally, note, that the electric constraints Y[ q Qq are taken care of by the product Y[ B p as 
had non-admissible labels met at a vertex, a corresponding q — 6j symbol gives zero. 
This concludes the proof of (fT5l l. 



For example we can introduce a black and white colour on boundaries and glue only black to white 
boundary, that is one with di weights to one with no weights of vertices. 
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